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We study a three-dimensional Anderson-Hubbard model under the coexistence of short-range 
interaction and diagonal disorder within the Hartree-Fock approximation. We show that the 
density of states at the Fermi energy is suppressed in the metallic phases near the metal- 
insulator transition as a proximity effect of the soft Hubbard gap in the insulating phases. The 
transition to the insulator is characterized by a vanishing density of states (DOS) in contrast to 
the formation of a quasiparticle peak at the Fermi energy obtained using the dynamical mean 
field theory in pure systems. Furthermore, we show that there exist frozen spin moments in the 
paramagnetic metal. 
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1. Introduction 

Understanding the nature of metal-insulator transi- 
tions (MIT) has been a central issue in condensed matter 
physics for a long time. 1 -' In particular, recently, the in- 
terplay of electron correlation and randomness has been 
attracting much attention experimentally and theoreti- 
cally because of their inevitable coexistence in real ma- 
terials. 

The single-particle DOS is a typical physical quantity 
that characterizes not only the nature of the MITs but 
also the electronic and magnetic properties in their vicin- 
ity. When the electron correlation causes the MIT as 
Mott transition, 2 ^ it opens a gap in the single-particle 
DOS. On the other hand, the Anderson insulator, which 
is driven by randomness, exhibits no gap in the single- 
particle DOS. 3 ' Such contrasting behavior of the DOS for 
these two types of MITs raise naturally a simple question: 
How does the DOS behave near Mott-Anderson transi- 
tions with coexisting electron correlation and random- 
ness? Despite extensive theoretical studies for several 
decades, 4 ' the nature of the Mott-Anderson transition 
has not yet been fully clarified. 

The Anderson-Hubbard model with coexisting on-site 
repulsion and diagonal disorder is one of the simplest 
models suitable for investigating the nature of the Mott- 
Anderson transition. Recently, we have determined the 
ground-state phase diagram of the three-dimensional 
Anderson-Hubbard model within the Hartree-Fock (HF) 
approximation. Furthermore, we have found an uncon- 
ventional soft gap (soft Hubbard gap) over the entire in- 
sulating phases. 5,6 ' Because only the short-range interac- 
tion is present in the Anderson-Hubbard model, the soft 
Hubbard gap cannot be explained by the conventional 
theory that attributes the formation of the soft gap to 
the long-range nature of the Coulomb interaction. 7 ' In- 
deed, we have proposed a multivalley energy landscape, 
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which may be characteristic of random systems, as the 
origin of the soft Hubbard gap. This observation of the 
soft gap is in clear contrast to the results of a numerical 
study within the dynamical mean-field theory (DMFT) 
8 ) and some mean-field studies 9, 10 ' that indicate the ab- 
sence of the soft gaps. This contradiction may be due to 
the insufficient treatment of spatial correlation, which is 
essential in the formation of the soft gap, in those studies. 

On the other hand, the behavior of the DOS in metals 
near MITs has been one of the central issues. In par- 
ticular, pseudogap phenomena observed in undcrdoped 
cupratc high-T c superconductors have inspired funda- 
mental discussions on the nature of the electronic states 
near the Mott insulator. n ) While the pseudogap and the 
soft gap both cause the reduction of the DOS at the 
Fermi level, the mechanism and origin of the pseudogap 
of the cuprates have not been established, and the role of 
randomness in stabilizing the superconductivity remains 
controversial. Recent numerical studies within the cellu- 
lar DMFT indicates the stabilization of a pseudogap 12 ) 
or Fermi arc 13-15 ' near the MITs in pure systems, in con- 
trast to the single-site DMFT results. 16 ' Because the soft 
gap mechanism may deepen and constructively stabilize 
the pseudogap formation in real experimental circum- 
stances with the inevitable coexistence of randomness 
and electron correlation, even on the HF level, further 
studies of the metallic phases near the Mott-Anderson 
transition will shed new light on the interplay of electron 
correlation and randomness and provide insight into the 
pseudogap phenomena. 

Spin polarization (or the formation of frozen spin 
moments) is another essential element in determining 
magnetic properties in the vicinity of the MITs, such 
as uniform magnetic susceptibility. Although a previous 
HF study on the three-dimensional Anderson-Hubbard 
model claimed the formation of frozen spin moments even 
in the paramagnetic metal, 17 ' the analyses are limited to 
finite system sizes and the bulk limit was not analyzed 
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after the extrapolation. 

In this paper, we show further numerical analyses of 
single-particle excitations and spin polarization on the 
metallic side within the HF approximation. 

2. Model and Method 

The Anderson-Hubbard Hamiltonian is defined as 

H = -t ]T c\ a c jc + uY, n n n n +2 W - ») n i° (!) 

on lattices with N s sites and N c electrons, where t is the 
hopping integral, U the on-site repulsion, c\ G (c, CT ) the 
creation (annihilation) operator for an electron with spin 
a on site i, rii C = c ia Ciai an d /i is the chemical potential. 
The random potential Vi is spatially uncorrelated and 
assumed to follow the Gaussian type distribution Py(Vi) 
with the average (V*) = 0: P v {Vi) = cxp(-lf /2cr 2 ) 

{a 2 = W 2 /12). We focus on half filling = U/2) on the 
cubic lattice throughout this paper. We take the lattice 
spacing as the length unit. 

We employ the Hartree-Fock (HF) approximation, 
where the wave function is approximated by a single 
Slater determinant consisting of a set of orthonormal 
single-particle orbitals {4> n } (n is an orbital index). The 
HF equation reads 

{Ho + Uy^X( n n} n n + (nit)nii)}<f>n = (2) 

i 

where is the one-body part of the Hamiltonian and 
we neglect (cLcij.) and (cj, c^). To find a site-dependent 
mean-field solution (rii a ) for the HF equations, we em- 
ploy the iterative scheme. One typically needs from sev- 
eral to several tens of initial guesses to obtain conver- 
gent physical quantities such as antifcrromagnctic (AF) 
order parameters and DOS. Here, wc employ pscudo-one- 
dimensional unit cells of L x L x M, where M 3> L. 

3. Results: Electronic and Magnetic Properties 
of Paramagnetic Metal 

3. 1 Density of states 

In Fig. 1, we present the calculated ground-state phase 
diagram within the HF approximation. 5,6 ' Hereafter, 
we take the hopping integral t as the energy unit. At 
U = 0, the Anderson- Hubbard model undergoes a metal- 
insulator transition (Anderson transition) from the para- 
magnetic metal (PM) to the paramagnetic insulator (PI) 
at a finite strength of disorder, W c = 21.29 ± 0.02. 18 ' 
On the other hand, at W = 0, since the system is half- 
filled with perfect nesting, the ground state is the anti- 
ferromagnetic insulator (AFI) for any nonzero value of 
U. Here, we discuss the ground-state phase diagram for 
U,W > 0. First, we focus on the spin degree of free- 
dom. For W > 0, the ground state is paramagnetic 
near U = 0. With increasing interaction, the ground 
states undergo an antiferromagnetic transition at a crit- 
ical point U c (> 0). Within the resolution of our calcula- 
tion, U c monotonically increases as disorder strength W 
increases. Next, we focus on the charge degrees of free- 
dom. The ground state is insulating for U, W 3> 1, which 
contains AFI as well as PI (PI is usually identified as an 
Anderson insulator) . Metallic phases identified by the di- 




Fig. 1. (Color online) Ground-state phase diagram of three- 
dimensional Anderson-Hubbard model at half filling for Gaus- 
sian distribution of P v . AFI, AF insulator; AFM, AF metal; PI, 
paramagnetic insulator (Anderson insulator); PM, paramagnetic 
metal. 



vergent localization length are restricted to a dome-like 
region (U < 6 and W < 25). 

Now, wc discuss the DOS in the paramagnetic metal. 
Figure 2 (a) shows the DOS at U = 2 (12 < W < 22) 
across the MIT point. At U = 0, the DOS shows no 
anomaly in the paramagnetic metal even near the MIT, 
because the DOS is not an order parameter of the An- 
derson transition and A(E-p) is nonzero in the Ander- 
son insulator. In contrast, for U > 0, the DOS at the 
Fermi energy should be zero in the paramagnetic in- 
sulator, because of the formation of the soft Hubbard 
gap. 5 ' 6 ' Indeed, the DOS shows a dip even in the metal- 
lic phase for 12 < W < 20 as a result of the proximity 
effect of the soft Hubbard gap. The dip gradually be- 
comes deeper with increasing W, and finally turns into 
the soft gap at the MIT. We have confirmed that the 
DOS for L = 8 and L = 10 are not different within the 
error bars at W — 14, indicating that the reduction of 
the DOS remains in the bulk limit. This critical behav- 
ior of the DOS is opposite to the results of a previous 
numerical study within the framework of the dynami- 
cal mean-field theory (DMFT) in infinite spatial dimen- 
sions, which support the divergence of A(Ep) toward the 
MIT. 8 ' This suppression of the DOS in our HF study 
indicates that the physical properties of the paramag- 
netic metal for U > is different from those of the non- 
interacting paramagnetic metal. For example, the Fermi 
surface in the paramagnetic metal vanishes as the dip be- 
comes deeper toward the MIT. Because the momentum- 
resolved DOS can be observed directly by angle-resolved 
photocmission spectroscopy, how the Fermi surface van- 
ishes in the momentum space is an interesting problem 
left for future studies. Furthermore, because the mag- 
netic susceptibility is proportional to the DOS at the 
Fermi energy, the Pauli susceptibility should be strongly 
reduced toward the MIT. However, it should be noted 
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Fig. 2. (Color online) (a) DOS with system size 8 X 8 X 250 at 
U = 2 (W = 22, 20, 18, 16, 14, and 12). We employ Lorentz 
broadening with a broadening factor of 5 X 10~ 4 . (b) The EA 
order parameter at U = 2 as a function of inverse length 1/L. 
(c) W dependence of the EA order parameter extrapolated to 
the bulk limit at U = 2. 



that the clarification of contributions from the incoher- 
ent part is needed to reach the complete understanding 
of the magnetic properties. 



the EA order parameter to the bulk limit. In fact, wc 
adopt M = 1000,250,250, and 100 for L = 4,6,8, and 
10, respectively. We find that the EA order parameter 
is extrapolated to a nonzero value even in the thermo- 
dynamic limit of the paramagnetic metal region, indi- 
cating the existence of frozen moments arising from the 
coexistence of the electron correlation and randomness. 
As shown in Fig. 2(c), the EA order parameter extrapo- 
lated to the bulk limit exhibits no singularity at the MIT. 
This formation of frozen moments in the paramagnetic 
metal is consistent with the result of a previous study 
using an effective-field theory. 20 ) However, the present 
HF calculations break the SU(2) symmetry and ignore 
quantum fluctuations. We need further studies beyond 
the mean-field level toward a complete understanding of 
the magnetic properties in the paramagnetic metal. 

4. Summary and Future Perspective 

We analyzed the three-dimensional Anderson- 
Hubbard model within the HF approximation. Wc found 
a dip in the DOS centered at the Fermi energy in metal- 
lic phases near the MIT for U > 0. The transition from 
a metal to an Anderson-Mott insulator is characterized 
by a continuously vanishing DOS in metals at the Fermi 
level, in sharp contrast to the picture of the dynamical 
mean field theory. Furthermore, we observed the forma- 
tion of frozen moments in the paramagnetic metal with 
coexisting electron correlation and randomness, which 
may be essential in determining the magnetic properties. 

In the previous papers, 5,6 ) as the origin of the soft 
Hubbard gap, we proposed the existence of low-energy 
multiply-excited states with electronic structures that 
arc globally relaxed from those of the ground state. There 
exist many excited states nearly degenerated with the 
ground state generically in a spin-glass phase. Therefore 
it will be interesting to clarify the connection between the 
spin-glass freezing and the formation of the soft Hubbard 
gap in the Anderson-Hubbard model. 

Before closing this paper, we make a brief remark con- 
cerning future topics. The present HF results support the 
stabilization of a pscudogap in a metal by the interplay 
of electron correlation and randomness, which may play 
a certain role in the formation of the pseudogap in vari- 
ous materials such as cuprates 21 * 1 and Cai-^Sr^VOa. 22 ) 
In the context of the Fermi arc observed in photoemis- 
sion for the cuprates, 21 ) it will be of great interest to 
investigate momentum-resolved DOS in a metal using a 
two-dimensional Anderson-Hubbard model. 
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3.2 Spin polarization/ formation of frozen spin mo- 
ments 

Next, we discuss spin degrees of freedom contained 
also in the incoherent part. The Edwards- Anderson (EA) 
order parameter for the spin glass is given by q = 
4A" 1 J2i (Si) 2 - 19 ^ Figure 2(b) shows an extrapolation of 



1) For a review see, M. Imada, A. Fujimori and Y. Tokura: Rev. 
Mod. Phys. 70 (1998) 1039. 

2) N. F. Mott: Proc. Phys. Soc. London Sect. A 62 (1949) 416. 

3) P. W. Anderson: Phys. Rev. 109 (1958) 1492. 

4) D. Belitz and T. R. Kirkpatrick: Rev. Mod. Phys. 66 (1994) 
261. 



4 J. Phys. Soc. Jpn. 



Letter 



H. Shinaoka et al. 



5) H. Shinaoka and M. Imada: Phys. Rev. Lett. 102 (2009) 016404. 

6) H. Shinaoka and M. Imada: J. Phys. Soc. Jpn. 78 (2009) 094708. 

7) A. L. Efros and B. I. Shklovskii: J. Phys. C 8 (1975) L49. 

8) V. Dobrosavljevic and G. Kotliar: Phys. Rev. Lett. 78 (1997) 
3943. 

9) V. Dobrosavljevic, A. A. Pastor, and K. Nikolic: Europhys. Lett. 
62 (2003) 76. 

10) K. Byczuk, W. Hofstetter, and D. Vollhardt: Phys. Rev. Lett. 
94 (2005) 056404. 

11) T. Timusk and B. Statt: Rep. Prog. Phys. 62 (1999) 61. 

12) Y. Z. Zhang and M. Imada: Phys. Rev. B 76 (2007) 045108. 

13) T. D. Stanescu and G. Kotliar: Phys. Rev. B 74 (2006) 125110. 

14) T. D. Stanescu, M. Civelli, K. Haulea, and G. Kotliar: Ann. 
Phys. (N.Y.) 321 (2006) 1682. 

15) S. Sakai, Y. Motome, and M. Imada: Phys. Rev. Lett. 102 
(2009) 056404. 



16) A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg: Rev. 
Mod. Phys. 68 (1996) 13. 

17) M. A. Tusch and D. E. Logan: Phys. Rev. B 48 (1993) 14843. 

18) K. Slcvin and T. Ohtsuki: Phys. Rev. Lett. 82 (1999) 382. 

19) S. F. Edwards and P. W. Anderson: J. Phys. F: Met. Phys. 5 
(1975) 965. 

20) M. Milovanovic, S. Sachdev, and R. N. Bhatt: Phys. Rev. Lett. 
63 (1989) 82. 

21) M. R. Norman, H. Ding, M. Randeria, J. C. Campuzano, T. 
Yokoya, T. Takeuchi, T. Takahashi, T. Mochiku, K. Kadowaki, 
P. Guptasarma, and D. G. Hinks: Nature 392 (1998) 157. 

22) R. Eguchi, T. Kiss, S. Tsuda, T. Shimojima, T. Mizokami, T. 
Yokoya, A. Chainani, S. Shin, I. H. Inouc, T. Togashi, S. Watan- 
abe, C. Q. Zhang, C. T. Chen, M. Arita, K. Shimada, H. Na- 
matame, and M. Taniguchi: Phys. Rev. Lett. 96 (2006) 076402. 



